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Abstract 

The proposal of pilgrim dark energy is based on speculation that 
phantom-like dark energy (with strong enough resistive force) can 
prevent black hole formation in the universe. We explore this phe¬ 
nomenon in loop quantum cosmology framework by taking Hubble 
horizon as an infra-red cutoff in pilgrim dark energy. We evaluate 
the cosmological parameters such as Hubble, equation of state pa¬ 
rameter, squared speed of sound and also cosmological planes like 
uj& — uo'fl and r — s on the basis of pilgrim dark energy parameter (u) 
and interacting parameter (d 2 ). It is found that values of Hubble pa¬ 
rameter lies in the range 744ooo5- It is mentioned here that equation 
state parameter lies within the ranges — 1 T 0.00005 for u = 2, 1 and 
(—1.12,—1), (-5,-1) for u = — 1,— 2, respectively. Also, 
planes provide ACDM limit, freezing and thawing regions for all cases 
of u. It is also interesting to mention here that — uA planes lie in 
the range (w# = — 1.131q2 5 , A? < 1-32). In addition, r — s planes 
also corresponds to ACDM for all cases of u. Finally, it is remarked 
that all the above constraints of cosmological parameters shows con¬ 
sistency with different observational data like Planck, WP, BAO, Hq 
and SNLS. 
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1 Introduction 

The accelerated expansion of the universe is one of the biggest achievements 
in the subject of cosmology [Tj. This expansion phenomenon follows through 
mysterious form of force called dark energy (DE). However, the nature of DE 
is still unknown. Different researchers have tried to explore the nature of DE 
through various aspects via theoretical and observational ways. As a result, 
they proposed different dynamical DE models as well as modified theories 
of gravity. The dynamical DE models have been developed in the scenarios 
of general relativity and quantum gravity. The pioneer candidate of DE is 
cosmological constant but it has two severe problems [2j. As an alternative 
to this candidate, the proposals of family of chaplygin gas [3], holographic 
BUS], new agegraphic [6] , polytropic gas |7j , pilgrim 0-ra DE models have 
been come forward. 

The holographic DE (HDE) has become an attractive DE model nowa¬ 
days, which is developed in the context of quantum gravity and widely used 
in solving the cosmological problems. The main idea of this model has come 
from holographic principle which is stated as the number of degrees of free¬ 
dom of a physical system should scale with its bounding area rather than its 
volume HU- With the help of this principle, a relationship between ultravi¬ 
olet and infrared (IR) cutoffs has been proposed by suggesting that the size 
of a system should not exceed the mass of black hole (BH) of the same size 
[12] , By using this relationship, Li [5] developed HDE density as follows 

p A = 3n 2 M 2 L~ 2 , 

here, n , M p . L indicate HDE constant, the reduced Planck constant, IR 
cutoff, respectively. On the basis of compatibility of HDE with the present 
day observations, different IR cutoffs have been proposed which includes 
Hubble, particle, event horizons, conformal age of the universe, Ricci scalar, 
Granda-Oliveros and higher derivative of Hubble parameter [13]- [15] etc. 

According to Cohen et al. [12], the bound of energy density from the idea 
of formation of BH in quantum gravity. However, it is suggested formation 
of BH can be avoided through appropriate repulsive force which resists the 
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matter collapse phenomenon. This force can only provide phantom DE in 
spite of other phases of DE like vacuum and quintessence DE. By keeping in 
mind this phenomenon, Wei |8]j has suggested the DE model called pilgrim 
DE (PDE) on the speculation that phantom DE possesses the large negative 
pressure as compared to the quintessence DE which helps in violating the null 
energy condition and possibly prevent the formation of BH. In the past, many 
applications of phantom DE exist in the literature. For instance, phantom 
DE is also play an important role in the wormhole physics where the event 
horizon can be avoided due to its presence [I6j. 

Also, it plays role in the reduction of mass due to its accretion process 
onto BH. Many works have been done in this support through a family of 
chaplygin gas H7j. It was also argued in the context of scalar field that BH 
area reduces up to 50 percent through phantom scalar field accretion onto 
it [18]. According to Sun [19], mass of BH tends to zero when the universe 
approaches to big rip singularity. It was also suggested that BHs might not 
be exist in the universe in the presence of quintessence-like DE which violates 
only strong energy condition [2D] - However, these works do not correspond 
to reality because quintessence DE does not contain enough resistive force to 
in order to avoid the formation of BH. 

The above discussion is motivated to Wei [8] in developing the PDE 
model. He analyzed this model with Hubble horizon through different theo¬ 
retical as well as observational aspects. Also, Saridakis et al. pTJ-JSUj have 
discussed the widely the crossing of phantom divide line, quintom as well 
as phantom-like nature of the universe in different frameworks and found 
interesting results in this respect. Recently, we have investigated this model 
by taking different IR cutoffs in flat as well as non-flat FRW universe with 
different cosmological parameters as well as cosmological planes [21 [ID]J . This 
model has also been investigated in different modified gravities [31|-|33j. In 
the present paper, we check the role of PDE in loop quantum cosmology 
(LQC). We develop different cosmological parameters and planes. The for¬ 
mat of the paper is as follows. In the next section, we provide the basic 
equations corresponding to PDE models. Also, we discuss the Hubble pa¬ 
rameter, EoS parameter and squared speed of sound in section 3. Section 
4 explores u# — oj# as well as statehnders planes. In the last section, we 
summarize our results. 
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2 Loop Quantum Cosmology and Pilgrim Dark 
Energy 


Nowadays, the discussion of DE phenomenon has also been done widely in 
the context of LQC to describe the quantum effects on our universe. The 
LQC is an interesting and attractive application of the Loop Quantum Grav¬ 
ity in the cosmological framework and it possesses the properties of non- 
perturbative and background independent quantization of gravity [33]-[39]. 
In recent years, many DE models have been studied in the scenario of LQC 
[30 , 131] . Jamil et al. [32] have explored cosmic coincidence problem phe¬ 
nomenon of modern cosmology by taking modified chaplygin gas coupled to 
dark matter. Also, some authors found that the future singularity appearing 
in the standard FRW cosmology can be avoided by loop quantum effects [33] . 
Chakraborty et al. |33J have made observational study of modified chaplygin 
gas in LQC. Here, we develop basic scenario of interacting PDE (with Hubble 
horzion) with cold dark matter (CDM) in LQC. The equation of motion in 
LQC has the form 


H 2 




( 1 ) 


where p indicates the sum of CDM and PDE densities. Also, p c = 167r 2 ^ 3 G 2 ^ 
represents the critical loop quantum density, 7 appears as dimensionless 
Barbero-Immirzi parameter. It is predicted that the big bang, big rip and 
other future singularities at semi classical regime can be avoided in LQC. 
Moreover, the modification in standard FRW cosmology due to LQC be¬ 
comes more dominant and the universe begins to bounce and then oscillate 
forever. 

It is argued that phantom DE with strong negative pressure can push 
the universe towards the big rip singularity where all the physical objects 
lose the gravitational bounds and finally dispersed. The PDE model is also 
developed in the favor of this scenario which is defined as 


p» = 3n 2 m A ~ u L~ u , (2) 

here u represents the PDE parameter. Wei explored the PDE model with 
different possible theoretical and observational ways to make the BH free 
phantom universe with Hubble horizon (L = // 1 ) through PDE parameter. 

In this work, we also choose PDE with the Hubble horizon which is the pi¬ 
oneer IR cutoff. Initially, it is plagued with a problem that its EoS parameter 
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provides inconsistent behavior with present status of the universe [IB]. This 
deficiency has been settled down with the passage of time by pointing out 
that HDE with this 1R cutoff can explain the present scenario of the universe 
in the presence of interaction with DM [45J. Also, the results of different cos¬ 
mological parameters have been established through different observational 
schemes by choosing HDE model with Hubble scale [IS mj. Sheykhi [1H] 
has discussed this model by taking interaction with CDM and pointed out 
that such model possesses the ability to explain the present scenario of the 
universe. 

In this work, we take interaction between PDE with CDM which takes 
the following form 

Pm + 3 H p m = 0, p$ + 3H(p# + p$) = —0, (3) 

where 0 possesses dynamical nature and appears as interaction term between 
CDM and PDE. Different forms of this interaction term has been proposed 
out of which we use the following form 

0 = 3 d 2 Hp m , (4) 

where d 2 is an interacting constant which appears as interaction parameter 
and exchanges the energy between CDM and DE components. This form 
of interaction term has been explored for energy transfer through different 
cosmological constraints. The sign of coupling constant decides the decay of 
energies either DE decays into CDM (when the interacting parameter is pos¬ 
itive) or CDM decays into DE (when the interacting parameter is negative). 
The present analysis from different aspects imply that the phenomenon of 
DE decays into CDM which is more acceptable and favors the observational 
data. Hence, the Eqs. (]3]) and (jl]) give 

Pm = p m0 a 3{d2 ~ 1) . (5) 

Also, by taking the differentiation of p$ (with Hubble horizon) with respect 
to x = In a, we get 

, H 

Pi) = U Pr>Jp- ( 6 ) 

3 Cosmological Parameters in LQC 

In this section, we will discuss the physical significance of cosmological pa¬ 
rameters corresponding to PDE with Hubble horizon in LQC scenario. 


5 


PDE with u=2 



Figure 1: Plot of H versus a for PDE in LQC with u = 2. 


PDE with u-1 



Figure 2: Plot of H versus a for PDE in LQC with u — 1. 

3.1 Hubble Parameter 

In order to check the behavior of Hubble parameter in this framework, we 
can find the following expression after some calculations 

H _ Pm 

IP ~ ~ MpH 2 (a) 

We solve the above differential equation (J7J) numerically by using Eqs.dTj)- 
(171) in terms of H and plot it against scale factor a for four different values 
of u = 2, 1, —1, —2 as shown in Figures 1 - 4 . The initial condition of H 
is taken as H[ao] — 74 as mentioned in Planck observations [49]. It has 
been greatly improved the precision of the cosmic distance scale through two 
recent analysis. Riess et ah [50] use HST observations of Cepheid variables 
in the host galaxies of eight SNe la to calibrate the supernova magnitude- 
redshift relation. Their best estimate of the Hubble constant, from fitting 


2 p c 


up$ 


p c - 2(/),? + p m ) 3 M 2 H 2 (a) 


(7) 
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PDE with u=—1 



Figure 3: Plot of H versus a for PDE in LQC with u = —1. 
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Figure 4: Plot of H versus a in LQC with u = —2. 
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the calibrated SNe magnitude-redshift relation, is 


H 0 = 73.8 ± 2Akms 1 Mpc 1 (Cepheids+SNe la), 


where the error is la level and includes known sources of systematic errors. 
Freedman et al. [51], as part of the Carnegie Hubble Program, use Spitzer 
Space Telescope mid-infrared observations to recalibrate secondary distance 
methods used in the HST Key Project. These authors find 

H 0 = [74.3 ± 1.5(statistical) ± 2.1 (systematic)] kms ~ 1 MpcT 1 

(Carnegie HP). 

In the present paper, it can be observed that the trajectories of Hubble 
parameter H(a ) attains the values approximately to 74+QQQg. Hence, the 
present results of Hubble parameter shows consistency with the above results 
found through observations. 


3.2 The Equation of State Parameter 

By using all above equations, we can obtain the EoS parameter as 


i j2 Pm . 
(jJ$ — —1 — (1 -h 


'lipr 


2 p c 


up# 


3 H 2 (a) [ p c ~ 2 (p& + p m ) 3M 2 H 2 (a) 


( 8 ) 


The plots of EoS parameter versus a are shown in Figures 5-8 for four dif¬ 
ferent values of u. In Figure 5 (u — 2), the trajectory of uj$ starts from 
phantom region and with the passage of time, it approaches to ACDM limit 
for the interacting cases d 2 = 0.2, 0.3. However, it remains in the phantom 
region for the interacting case d 2 = 0.4. In case of ii — 1 (Figure 6), the tra¬ 
jectories of EoS parameter remains in quintessence region for d 2 = 0.4 while 
it approaches (from quintessence region) to ACDM limit for the other two 
cases of d 2 . For u — — 1, —2 (Figures 7-8), the EoS starts from phantom with 
comparatively high value and goes towards ACDM limit for d 2 = 0.2, 0.3 and 
always remains in phantom for other case of d 2 . Moreover, the constraints 
on EoS parameter has been put forward by Ade et al. |49] (Planck data) 







PDE with u=2 



Figure 5: Plot of versus a for PDE in LQC with u = 2. 


PDE with u= 1 



Figure 6: Plot of co# versus a for PDE in LQC with u = 1. 
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Figure 7: Plot of u># versus a for PDE in LQC with u = —1. 
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PDE with u=-2 
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Figure 8: Plot of versus a for PDE in LQC with u = —2. 


LOfl = 

1 1 Q+0.24 
■ L ' io -0.25 

(Planck±WP±BAO), 

UJ$ = 

-1.09 ±0.17, 

(Planck±WP±Union 2.1) 

= 

1 1 Q+0.13 
± - io -0.14> 

(Planck±WP±SNLS), 

UJ$ = 

—1 24 +0 ' 18 

(Planck±WP±tf 0 ), 


by implying different combination of observational schemes at 95% confidence 
level. It can be seen from Figures 5-8 that the EoS parameter also meets the 
above mentioned values for all cases of interacting parameter which shows 
consistency of our results. The above discussion shows that all the models 
provides fully support the PDE phenomenon. 


3.3 The Square Speed of Sound 

In order to analyze the stability of PDE model in this scenario, we extract 
the squared speed of sound which is given by 


2 P P 

v s = - = - r 

P P 


(9) 
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PDE with u=2 



Figure 9: Plot of v 2 versus a for PDE in LQC with u = 2. 


where pressure corresponds to PDE only. After some calculations, we can 
obtain squared speed of sound as follows 

v 2 = i(—3 - (a~ 3+3d2 d 2 p m0 H(a)~ u )n~ 2 + (a~ 3+3d2 p m0 u(-un 2 

x H(a)~ 2+U + 2p c (-2a~ 3+3d2 d + p c - 6 n 2 H{a) u )~ 1 )) 
x H~ 2 (a) — (a~ 3+3d2 d 2 p m0 u 2 — (3( —1 + d 2 )d 2 H (a) 4 - “)n ~ 2 

- 3(—1 + d 2 )u 2 n 2 H(a) 1+u - a~ 3+3d2 p m 0 (-3 + u)u 3 n*H(a)- 3+2u 
+ ( 6 a 6 (—1 + bs)up c H(a) 3 (p c — 6n 2 H(a) u ))(2a 3d ~ p m0 — a 3 a 

+ 6 a 3 n 2 H(a) u )~ 2 + (2a 3d2 d 2 p m oup c H(a) 2 ~ u )(n 2 (2a 3d2 p m0 — a 3 p c 
+ 6 a 3 n 2 H{a) u ))~ l + (2a 3d2 p m0 u 2 p c n 2 H (a)~ 1+ “(—(—4 + u)(2a 3d2 p m0 

— a 3 p c ) + 24a 3 n 2 H(a) u ))(2a 3d2 p m0 — a 3 p c + 6 a 3 n 2 H (a) u )~ 2 

+ (4a 3+3d2 p m0 up 2 c H(a)(2a 3d2 p. m0 - a 3 p c + 6 a 3 (l + u)n 2 H(a) u ))(2a 3d2 
x pmo — a 3 p c + 6 a 3 n 2 H (a) u )~ 3 )(ciuH (a) 3 (—iin 2 H (a)~ 2+u 
+ (2p c )(-2a~ 3+3d2 p m0 + p c - 6n 2 H(a) u )~ 1 ))~ 1 ). 

The plots of squared speed of sound versus a for three different values of 
d 2 and four values of u — 2, 1, —1, —2 is shown in Figures 9 - 12 . ft can be 
observed from Figures 9-10 (for cases u = 2, 1) that the squared speed of 
sound remains negative for all cases of d 2 which exhibits the instability of the 
PDE in LQC scenario. For the cases u = —1, —2 (Figures 11 - 12 ), it exhibits 
the stability of the present model for all cases of d 2 . 
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PDE with u-1 



Figure 10: Plot of v 2 versus a for PDE in LQC with u — 1. 


PDE with u=— 1 



Figure 11: Plot of v 2 versus a for PDE in LQC with u = —1. 


PDE with u=-2 



Figure 12: Plot of v 2 versus a for PDE in LQC with u = —2. 
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4 Cosmological Planes in LQC 

Here, we will discuss the physical significance of cosmological planes corre¬ 
sponding to PDE with Hubble horizon in LQC scenario. 

4.1 (jjtf — Analysis 

Here, we find the regions on the u>$ — co# plane ( co# represents the evolution of 
u)&) as defined by Caldwell and Linder [52] for models under consideration. 
The models can be categorized in two different classes as thawing and freez¬ 
ing regions on the — co$ plane. The thawing models describe the region 
co# > 0 when co# < 0 and freezing models represent the region co# < 0 when 
Utf < 0. Initially, this phenomenon was applied for analyzing the behavior of 
quintessence model and found that the corresponding area occupied on the 
— u$ plane describes the thawing and freezing regions. Differentiating co# 
with respect to x — In a and after some calculations, we obtain 

ui# = (3iL(o) 4 ) _ 1 o _ 4 + 3 d 2 p m 0 (o _ 3 + 3 d 2 (i 4 p m on 2 — (3(—1 + d 2 )d 4 H(a ) 4: ~ u ) 
x nr 2 — 3( —1 + d 2 )u 2 n 2 H (a) 1+u — a~ 3+3d2 p m0 ( — 3 + u)u 3 n A H(a)~ 3+2u 
+ (6a 6 ( —1 + d 2 )up c H(a) 3 (p c — 6n 2 H (a) u ))(2a 3d ~ p m0 — a 3 p c + 6a 3 
x n 2 H(a) u )~ 2 + (2a 3d2 d 4 p m0 up c H (a) 2 ~ u ){n 2 (2a 3d “p m0 - a 3 p c 
+ 6 a 3 n 2 H(a) u )Y 2 + (2a 3d2 p m ou 2 p c n 2 H(a)~ 1+u (-(-4: + u)(2a 3d2 p m0 

— a 3 p c ) + 24a 3 n 2 H(a) u ))(2a 3d2 p m0 — a 3 p c + 6 a 3 n 2 H (a) u )~ 2 + (4a 3+3rf2 
x p m0 up 2 c H(a)(2a 3d2 p m0 - a 3 p c + 6 a 3 (l + u)n 2 H(a) u ))(2a 3d2 p m0 

— a 3 p c +Ga 3 n 2 H(a) u )~ 3 ). 


We also construct the u$ — co 1 # plane for PDE model with different values 
of u in LQC as shown in Figures 13 - 16 . In all cases of u and d 2 , the uj$ — co# 
plane corresponds to ACDM limit, i.e., (cotf,co#) = (—1,0). Also, u& — co# 
plane shows thawing regions for the cases u = 2, —1, —2 (Figures 13 , 15 , 16 ) 
and corresponds to freezing region for the case u — 1 (Figure 14 ). It has 
been developed the following constraints on w# and iu# by Ade et al. [39]: 

, _ _11 q+0.24 

U)fl — I.Io_q25 

co< 1.32, 


(Planck+WP+BAO), 

(Planck+WP+BAO) 
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PDE with u=2 



Figure 13: Plot of — u >^ for PDE in LQC with u = 2. 


PDE with u=l 



Figure 14: Plot of u# — for PDE in LQC with u = 1. 


PDE with u=—1 



Figure 15: plot of oj\ — u' A for PDE in LQC with u — — 1. 
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PDE with u=—2 



Figure 16: plot of — oj' a for PDE in LQC with u = —2. 


at 95% confidence level. Also, other data with different combinations of ob¬ 
servational schemes such as (Planck+WP+Union 2.1) and (Planck+WP+SNLS) 
favor the above constraints. In the present case, the trajectories of cj# against 
ujtf also meet the above mentioned values for all cases of interacting parame¬ 
ter which shows consistency of our results as shown in Figures 13 - 16 . Hence, 

— ui'tf plane provides consistent behavior with the present day observations 
in all cases of u. 


4.2 Statefinder Parameters 

The statefinder parameters are defined as follows [53] 

a r — 1 

r = TUP’ S= 3 {q-\y 


( 10 ) 


where q is the deceleration parameter. These parameters are dimensionless 
and possess the ability to explain the current accelerated scenario. These 
parameters have geometrical diagnostic due to their total dependence on the 
expansion factor. The statefinders are useful in the sense that we can find the 
distance of a given DE model from ACDM limit. The well-known regions 
described by these cosmological parameters are as follows: (r, s) = (1,0) 
indicates ACDM limit, (r, s ) = (1,1) shows CDM limit, while s > 0 and r < 1 
represent the region of phantom and quintessence DE eras. By following the 
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PDE with u=2 



Figure 17: Plot of r — s for for PDE in LQC with u — 2. 
papers PHD], we can obtained the following form of statefinders 

r = 1 + ^a _3+3<i2 p m0 (— u 2 n 2 Ha~ 3+u — (d A H (a)~ u )n~ 2 + (2up c )(H(a) 

x (~2a~ 3+3d2 pmo + Pc~ < on 2 H(a) u ))~ l )(-3 - (a~ 3+3d ~ d 4 p m0 H(a)~ u ) 
x n~ 2 + (a~ 3+3d2 p m0 u(-un 2 H(a) ~ 2+u + (2p c )(-2a~ 3+3d2 p m0 + p c 
- Qn 2 H(a) u )~ 1 ))H(a)~ 1 - (3 a 3 ~ 3d2 H(a)oj' i) ){p m0 u{-un 2 H{a)~ 2+U 
+ (2p c )(-2a~ 3+3d2 p m0 + p c ~ 6n 2 H (a) u )~ 1 ))~ 1 ). (11) 

and 

s = -a~ 3+3d2 p m Q(—u 2 n 2 H(a)~ 3+U — (d i H(a)~ u )n~ 2 -\- {2up c )(H(a) 

9 

x (-2 :a - 3 + 3d2 p m0 + p c - 6n 2 H(a) u ))~ 1 )(-3 - (a' 3+3d2 d 4 p m0 ^(a) _ “) 
x n -2 + (a _3+3d2 p m0 n(— un 2 H(a)~ 2+w + (2p c )(—2a*- — 3 + 3 d 2 )p m0 
+ Pc~ Qn 2 H(a) u )- 1 ))H(a)- 1 - (3a 3 ~ 3d2 H{a)^)(p m0 u(-un 2 H{a)- 2+u 
+ (2p c )(-2a~ 3+3d2 p m0 + p c - 67i 2 H(a) u )~ 1 ))- 1 ). (12) 

We also develop r — s planes corresponding to the present cosmological 
scenario for different values of u as shown in Figures 17-20. The r — s 
corresponds to ACDM limit for all cases of u. 
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PDE with u= 1 



Figure 18: r — s for PDE for PDE in LQC with u — 1. 


PDE with u=— 1 



Figure 19: r — s for PDE in LQC with u = —1. 


PDE with u=-2 



Figure 20: r — s for PDE in LQC with u = —2. 
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5 Concluding Remarks 

We have considered the framework of interacting PDE with Hubble horizon 
in LQC framework. The main motivation of this work is to analyze the 
cosmological scenario as well as construction of possible constraints of PDE 
parameter u where it fulfill the PDE phenomenon. For this purpose, we 
have constructed Hubble parameter, EoS parameter, squared speed of sound, 
and r — s planes, numerically. We have discussed these parameters 
corresponding to four values of u — 2 , 1 , — 1 , —2 and three value of d 2 = 
0.2, 0.3, 0.4. We have observed that the trajectories of Hubble parameter 
H(a ) for all cases of u attains the values approximately to 74zto;8o5 (Figures 
1-4). These obtained range of H (a) shows consistency with the observational 
values such as H 0 = 73.8 ± 2.4 [50] and H 0 = 74.3 ± 1.5 [5Tj . 

Moreover, the EoS parameter also shows consistency with the present 
day observations. For instance, the trajectory of u;# exhibits the ranges 
—1 — 0.0050 and —1 + 0.00005 for the cases u — 2, 1 as shown in Figures 
5-6. For u = — 1, — 2 (Figures 7-8), the EoS parameter lies in the ranges 
( — 1.12, —1) and (—5, —1), respectively. These constraints on EoS parameter 
compatible with the constraints as obtained by Ade et al. [39] (Planck data) 
which is given as follows: 


= 

1 1 Q+0.24 
± - ±o -0.25 

(Planck±WP±BAO), 

U)ti = 

-1.09 ±0.17, 

(Planck±WP±Union 2.1) 

Ld$ = 

1 1 Q+0.13 

(Planck±WP±SNLS), 

UJ$ = 

—1 24 +0 ' 18 

(Planck±WP±77 0 )- 


The above constraints has been obtained by implying different combination 
of observational schemes at 95% confidence level. 

It can also be observed from Figures 9-10 (for cases u — 2, 1) that the 
squared speed of sound remains negative for all cases of d 2 which exhibits the 
instability of the PDE in LQC scenario. For the cases u = —1, —2 (Figures 
11 - 12 ), it exhibits the stability of the present model for all cases of d 2 . We 
have also observed that the c — oj# plane corresponds to ACDM limit, i.e., 
(u$, U#) = (—1, 0) in all cases of u and d 2 . Also, — plane shows thawing 
regions for the cases u — 2, —1, —2 (Figures 13 , 15 , 16 ) and corresponds to 
freezing region for the case u = 1 (Figure 14 ). It has been developed the 
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following constraints on w& and w ^ by Ade et al. [45]: 


= —1.1 3^25 (Planck+WP+BAO), 

u# < 1.32, (Planck+WP+BAO) 

at 95% confidence level. Also, other data with different combinations of ob¬ 
servational schemes such as (Planck+WP+Union 2.1) and (Planck+WP+SNLS) 
favor the above constraints. In the present case, the trajectories of u# against 
calso meet the above mentioned values for all cases of interacting parame¬ 
ter which shows consistency of our results as shown in Figures 13-16. Hence, 
c— l <4 plane provides consistent behavior with the present day observations 
in all cases of u. Also, the r — s corresponds to ACDM limit for all cases of u. 
Finally, it is remarked that all the cosmological parameters in the scenario 
of LQC with PDE shows compatibility with the current observations. 
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